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Abstract 

Fenimore et al. (1995) have recently presented a very tight correlation between the 
spectral and the temporal structure in Gamma Ray Bursts (hereafter GRBs). In particular, 
they discovered that the durations of the constituent subpulses which make up the time 
profile of a given GRB have a well defined power-law dependence, of index ~ 0.45, on the 
energy E of the observed photons. In this note we present two simple models which can 
account in a straightforward fashion for the observed correlation. These models involve: 
(a) The impulsive injection of a population of relativistic electrons and their subsequent 
cooling by synchrotron radiation, (b) The impulsive injection of monoenergetic high energy 
photons in a medium of Thompson depth tt ~ 5 and their subsequent downgrading in 
energy due to electron scattering. We present arguments for distinguishing between these 
two models from the existing data. 
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I. Introduction 



The issue of the nature of GRBs, a contentious one since their discovery, has only 
deepened with the more recent observations by BATSE which confirmed their isotropic 
distribution in the sky while indicating a significant lack of faint events, inconsistent with 
a uniform distribution in space (Meegan et al 1992). These facts have provided support 
to models invoking the cosmological origin of the GRBs, thus increasing further the un- 
certainty associated with the location, nature and origin of the GRB phenomenon. This 
uncertainty is due largely to the lack of association of any class of interesting astrophysical 
objects within the positional error boxes of GRBs, which would set the distance scale and 
hence the luminosity associated with these events. In the absence of such positional asso- 
ciation, attempts to gain understanding on the nature of GRBs have relied in searches for 
systematics in the time and/or energy domain of the observed bursts. Part of our contin- 
uing ignorance with respect to the nature of GRbs is due to the absence of any significant 
correlations between the GRB spectral and/or temporal properties (Paczynski 1992). 

This situation seems to have recently changed: Fenimore et al. (1995) have presented a 
very well defined correlation between the average duration of the subpulses usually present 
in the time profile of a given burst - or the duration of the entire burst in the absence 
of subpulses - and the energy channel of observation. They have found that the average 
pulse width, Ar, as determined either by the autocorrelation function of the entire GRB 
or by specific fits of the time profiles of individual subpulses in an given profile, is well 
fitted by a power law in the energy E at which the observation is made; more specifically 
they have found that Ar ~ ^4^;-o.45±o.05^ w ] iere A is a, proportionality constant. The 
above reference thus made quantitative an effect that had been noticed earlier by Fishman 
et al. (1992) and Link, Epstein & Priedhorsky (1993) who indicated that, in general, the 
various peaks in the GRB time profiles are shorter and better defined at higher energies. 
This unique piece of information may, therefore, serve as a stepping stone in probing the 
physical processes associated with the radiation emission in GRBs and maybe the entire 
GRB phenomenon. More recently, Norris et al. (1996) have further elaborated on the time 
profiles of GRBs as a function of the energy channel of observation. These authors have 
found that the relation reported by Fenimore et al. (1995) indeed holds but it may have 
slightly different slopes depending on the overall position of the observed pulse within the 
entire burst. 

In this note we present two generic models which can provide an account of the par- 
ticular functional form of the observed correlation between Ar and E. It is of interest to 
note that the physical processes associated with these two specific models are fundamen- 
tally different, despite the fact that they both can produce the observed pulse duration - 
energy correlation. However, these models are quite simple in their concept and specific 
enough to consider possible that further observations and detail fits would yield systemat- 
ics which could discriminate between them or will lead to further insights into the radiation 
mechanism of GRBs. The first model involves the impulsive injection of a power law of 
relativistic electrons in a given volume and their subsequent cooling by synchrotron or 
inverse Compton radiation, while the second one the impulsive injection of high energy 
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photons in an cold medium of Thompson depth tt ~ 10 and their subsequent downgrading 
by electron scattering. 

In section II the mathematical formulation of both models is presented along with 
sample time-dependent spectra and model pulse durations at various frequencies. In sec- 
tion III the results are discussed and conclusions are drawn concerning the impact of these 
results on GRBs models. 

II. The Models 

To begin with, it is worth noting that the relation discussed by Fenimore et al. (1995) 
precludes from the outset the possibility that the observed spectra are the result of up- 
Comptonization of soft photons by repeated scatterings in a steady-state, hot electron 
population (see e.g. Sunyaev & Titarchuk 1980); in this case the the higher energy photons 
are those which would have spent longer time within the hot, scattering medium, leading 
to a correlation between the burst duration At and the photon energy E of opposite sense 
to that observed. 

Before we discuss the models for the energy dependence of GRB duration in detail, we 
present some simple, heuristic arguments which qualitatively indicate that a well defined 
spectro - temporal relation similar to that observed is not indeed at all surprising. These 
considerations are motivated by the recent cosmological scenario for GRBs suggested by 
Meszaros & Rees (1993) and also by some arguments concerning photon scattering in a 
thick media and should generally be as applicable to galactic GRB models. It is worth 
noting here that, despite their apparent dissimilarity, both models outlined below describe 
essentially the same process, namely the energy degradation of a population of relativistic 
particles, with one of the models refering to relativistic electrons while the other to photons. 
What is more important to bear in mind is that the observed correlation appears to require 
the presence of some sort of "cooling" in GRBs. 

Consider the impulsive injection of relativistic particles at a given volume and assume 
the emission to be optically thin. If synchrotron or inverse Compton losses dominate the 
evolution of the electron distribution, then the loss rate per electron is 7 oc 5 2 7 2 , where 7 is 
the Lorentz factor of a given electron and B the magnetic field (if the losses are dominated 
by inverse Compton, then B 2 should be replaced by the ambient photon energy density). 
Therefore, the characteristic electron life time is At ~ 7/7 oc 7 . Since the characteristic 
energy of synchrotron (or inverse Compton) emission is E oc 7 2 eo (^0 is a characteristic 
energy), expressed in terms of E rather than 7, the time duration is At oc E~ x l 2 . 

Consider, alternatively, the impulsive injection of monoenergetic high energy photons 
of energy Eq in a medium of Thomson depth tt 3> 1 and of electron temperature kT <C 
Eq. Photons of energy E suffer fractional energy loss of order AE/E = (1 — fj,)E/m e c 2 
per collision, where \i is the cosine of the photon scattering angle, or in terms of the 
photon wavelength A, AA = (1 — fi). Averaged over solid angle, (A A) = 1. So after n 
scatterings, (A n ) = n + Xq where Ao is the original wavelength (omitted for n > 1). In the 
diffusion regime, the dispersion of photon wavelengths is also propotional to the number 
of scatterings. Hence the following relation between the number of scatterings or time (in 
units of scattering time): At oc A 2 oc E~ 2 . For a finite medium, this diffusion process 



3 



lasts for a number of scatterings n oc r T 2 , thus the maximum width (in wavelegth) of 
the escaping photons would be AA oc tt- In this case, however, the precise law between 
the pulse width and the energy depends on the photon source distribution within the 
scattering medium, as it will be exhibited in more detail in the next section. Nonetheless, 
the qualitative correlation between the energy and duration of pulses will always be as 
that outlined above. 

In the following, we make the above heuristic arguments more concrete by presenting 
well formulated solutions of the problems outlined above. 

77 a. Relativistic Electron Cooling 

Consider the simplest case of time-dependent synchrotron radiation by a population 
of non-thermal electrons: Relativistic electrons of Lorentz factor 7, with a power law 
distribution of index T i.e. Q e (7) = Q7 _r \ are injected in a volume of tangled magnetic 
field B at an instant t = to and are left to cool by emission of synchrotron (on inverse 
Compton radiation). The equation which governs the evolution of the electron distribution 
function iV(7,t) in time is 

^ = -^mN(l,t)} + Q e ( 1 )8(t-t ) 
= ^-[(3^N{ 1 ,t)]+Q e { 1 )d{t-to) 

where (7) = — fi'j 2 is the mean electron loss rate and (3 = (4/3)(5 2 /87r)(oY/m e c 2 )c, 
while Q e {l)i as given above, represents the injection of relativistic electron population. 
This equation can be solved by the standard technique of converting it into an ordinary 
differential equation in t and integrating it along the characteristic curves 7 = ~f(t). Thus, 
multiplying through by 7 2 one obtains 

- /37 2 ;f [7 2 iV(7, t)] = 7 2 Qe(7)5(t - to) 
m 01 (la) 

^ = 7 2 Q e (7)^-t ) 

with the last step becoming obvious on identifying /dt with — P'j 2 . The total differential 
of equation (la) can be integrated along the integral curves of the electron energy 7; these 
are given by the solution of d'j/dt = —f3^ 2 , namely 7 = 70/(1+70/^), where 70 is the initial 
electron energy. Assuming Q e (l) = <11~ T \ the resulting form of the electron distribution 
function is 

AT/ry +) — f <?7~ r (l - I3^t) v ~ 2 for 7 <l//3t ( ) 

l7,tJ "\0 for 7 >l//ft • [Z) 
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Given the electron distribution function iV(7,t), the emitted spectrum can be calcu- 
lated by convolving AT (7, t) with the synchrotron emissivity 6(7, v) (see e.g. Rybicki & 
Lightman 1979). However, the form of the emerging spectrum can be easily calculated us- 
ing the 5-function approximation to the single electron emissivity, i.e. 6(7, v) oc 8(y — 7 2 ^ c ), 
where v c ~ 4x 10 6 B(G) Hz is the cyclotron frequency. Assuming, further, that the injected 
distribution function has a minimum Lorentz factor 7 m (i) the resulting spectrum has the 
form 



where v m (t) is the synchrotron frequency corresponding to the lowest energy electrons of 
the injected electron Lorentz factor 7 m (which is also a function of time), and v c is the 
cyclotron frequency corresponding to magnetic field B. 

Figure 1 shows the time dependent spectra obtained by integrating the distribution 
function of Eq. (2) over the exact expression of the synchrotron emissivity, with an injection 
spectrum oc 7~ 2 between Lorentz factors 7 m = 10 6 and 7m = 10 9 . These specific values 
(especially that of 7 m ) were chosen to assure that the low energy turnover of the spectrum 
occurs at v m ~ 10 20 Hz (= 414 keV), as indicated by the GRB spectra (Schaeffer et 
al. 1994), and that the duration of the burst be r& ~ a few seconds, also in agreement with 
observations. These considerations then serve to also determine the value of the magnetic 
field associated with this emission to be B ~ 30 Gauss. 

If the observed durations are indeed only apparent rather than intrinsic, the result of 
emission from a plasma in relativistic motion toward the observer, as suggested e.g. in the 
currently popular relativistic blast wave GRB models (see e.g. Meszaros & Rees 1992), 
then both the break frequency v m and the characteristic time scale At should be shifted by 
the appropriate powers of the bulk Lorentz factor of the flow Tto reflect the values at the 
rest frame of the fluid. In this case, the intrinsic frequency of emission, i/' m , will be shorter 
from the observed one z/ m o by a factor T, i.e. v' m ~ z/ m o /r while the intrinsic duration At' 
would be longer than the observed one At ~ a few sec, by a factor S~ x = T(l — (3l) — T, 
i.e., At' ~ At ■ V (f3 L is the velocity of the fluid, i.e. T = (1 — /3l) -1 ^ 2 )- Using the primed 
(rest frame) values of the duration and break frequency to determine our values of B and 
7 m leads to the following scalings for the magnetic field B and the value of 7 m on the rest 
frame of the fluid in terms of the observed duration Ar(sec) and z/ mj0 = 10 20 z/ 2 o Hz 




for v < v m (t) = 7 m (t) 2 z/ c 

for V > Vm(t) 
for v > vjf3 2 t 2 



(3) 



B ~ 10 5 / 3 r- 1 / 3 z/ 20 1/3 Ar- 2 / 3 (sec) Gauss 



and 



lm — 



10 17/3 r -l/3 z/ l/3 Ar l/3 (sec) 
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So while our calculations were performed with a field value of B ~ 30 Gauss and 
7 m ~ 10 6 the values of these parameters in the fluid rest frame would depend on the value 
of the bulk Lorentz factor V by the relations given above. 

Figure 2 shows the synchrotron emission as a function of time at a set of frequencies, 
spaced logarithmically, with each curve corresponding to a frequency x 10 smaller than the 
previous one. It is apparent in this figure that the burst duration is inversely proportional 
to the square root of the radiation frequency, in accordance with the relation reported by 
Fenimore et al. (1995). It is also apparent in this figure that, for the particular form of 
electron injection considered in this example (i.e. impulsive injection), there is no time lag 
in peak emission between frequencies in the range v > 7 m (£o) 2 ^c, i-e. for frequencies higher 
than the synchrotron frequency of the lowest energy electrons at the time of injection. Time 
lags in peak emission develop only for frequencies smaller than that above, resulting from 
the "cooling" of the lowest energy electrons to energies 7 m (£) < 7™(£o)- 



II b. Energy Downgrading by Electron Scattering 

Consider alternatively the injection of photons of energy E ~ m e c 2 near the center 
of a cold (T e <C E), spherical (for simplicity) electron cloud of Thomson depth tt ~ 5. 
As the photons random walk out of the cloud, they also lose energy in collisions with 
the ambient electrons. A monoenergetic photon pulse at the center of the cloud spreads, 
therefore, both in energy and in time. The duration of the pulse for photons of a given 
energy can be computed from the following considerations: 

At a given scattering event the fractional change of the photon energy is 



v m e c 

or, expressed in terms of the photon wavelength A = m e c 2 /hiy, 



Au = (1 - n)hv ^ 



AA=l-/x, (6) 

where fi is the cosine of the angle between the incident and scattered photon momenta. 

Let (f>(n) = (3/167r)(l + /U 2 ) denote the phase function of the Thompson cross section, 
such that J (p(p)dQ = 1; then the average change in the photon wavelength (AA) and the 
dispersion in the photon energies a\ per scattering are 

mx = (AA) = J 0(/x)(l - fi)dQ = 1 (7) 

and 

^ = ((AA) 2 )-(AA) 2 = ((AA) 2 >-1 

= j m^-^) 2 d^-i = lj\i+^)(i-^ 2 d^ = 2 - (8) 
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The distribution (probability) of photons over number of scatterings n with a mean 
number of scatterings u is given by a Poisson distribution, P u (n), which is the limit of 
binomial distribution for the case of rare events. For the Poisson distribution the values 
of the mean m n and dispersion cr^ are equal, i.e. a\ = m n = u. 

Now, the photon distribution in wavelength A for a given mean number of scatterings 
u (or alternatively after time u measured in units of the scattering time r = l/o~Tn e c) is 
given by the convolution 

oo 

J u (\) = J2*n(\)P u (n) (10) 

where ^ n (X) is the photon distribution after n scatterings, with mean wavelength (A n ) = 
n + Ao (Ao is the wavelength of the initial photon) and dispersion a x (n) = n- o\ = (2/5)n 
and P u (n) is the probability of n scatterings within time u. 

The mean of this distribution is 

/oo 00 / />oo \ 

J u (X)d\ = J2[ ^n(X)dX P u {n) 
•0 n= o / 

00 v / 

= ^2(n + X )P u (n) = u + A 

n=0 

and the dispersion 



/•OO 

°l{u)= / (X-m u x ) 2 J u (X)dX 

/•OO 

= / [(\- m lf + 2(X-m n x )(m n x -m u x ) + (m n x -m u x ) 2 ]dXP u (n) 



'A 

00 

E 

n=0 



-n + 0+ (n - uf 



(12) 



P u (n) 



2 7 
= —u + u = -u 
5 5 



Therefore the photon spectrum after time u will have the form 

(X-ml) 2 



MX) = 



1 



exp 



2a\{u) 



(13) 



which represents a gaussian of width u about A — Ao- In fact, the above relation can be 
viewed both as the photon spectrum escaping after time u or alternatively as the time 
history for photons of a given wavelength A. One can use the above expression to obtain 



7 



the dispersion in escaping time for photons of a given energy. This is simply given by the 
width of the gaussian in equation (13), which, considering the definitions of cr\(u) and m^, 
yields A — A ^ 2y/2u, i.e. that At oc E~ 2 , as suggested by our heuristic arguments. 

The emergent spectrum of the photons escaping over the duration of the burst is given 
by the convolution 

/•OO 

/(A) = / J u {\)P{u)du (14) 
Jo 

where P(u) is the probability of a photon to escape after u scatterings. As a concrete 
example, we consider a photon source distribution according to the first eigenfunction of 
the diffusion operator; with this assumption, the function P(u) has the form 



p(u) = p e -P« (15) 

with (3 = 7t 2 /(3(to + 2/3) 2 ) being the eigenvalue of the same operator. Then the integral 
of Eq. (14) above is analytic yielding 

/(A) = (1 + ^)1/4(1 + ^)1/10 ^^ - A )/3] (16) 

which indicates that the spectrum is a power law in energy of index —2 with an exponential 
turn-over at A ~ 1/(3. 

The analytic results of the above discussion should be viewed only as an illustrative 
example to indicate the relationship between the photon energy E and the duration of a 
photon pulse At. It is also apparent from the above discussion that the precise form of 
this relation depends on the distribution of photon sources within the scattering cloud, 
the relation derived from equation (13) correponding to a homogeneous, infinite medium. 
Our analytic calculations are also more appropriate for the case in which the high energy 
photon source at the center of the electron cloud. This ensures that the escaping photons 
will undergo quite a few scatterings before escape and that they will lose a significant 
fraction of their energy in the first few scatterings. These scatterings reduce their energy 
to the point where the treatment of radiative transfer can be replaced by diffusion in energy 
and space (i.e. small energy change per scattering at the Thomson cross section). Clearly, 
the random walk of photons situated close to the cloud boundary, cannot be described 
by diffusion with small energy change at each scattering event. A large fraction of the 
high energy photons (assuming monoenergetic high energy injection) undergoes only a few 
scatterings before escape but with large energy change per scattering (Hua & Titarchuk 
1996). Thus the accurate consideration of the photon transport distributed over a finite 
plasma cloud could, conceivably, lead to an energy - duration relation in closer agreement 
to the results of Fenimore et al. (1995). 

Because our analytic caclulations are appropriate only when the source of photons is 
located at the center of the cloud, we have extended our model to make it more realistic 
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at the expence of introducing an additional parameter. The extended model consists of 
two concentric spheres; it is assumed that cold electrons are uniformly distributed over 
the entire volume of the larger sphere of radius Rq. The smaller sphere, of radius R s , 
represents the volume within which the high energy photon sources are located and the 
photon injection takes place, uniformly over its entire volume. This model has an additional 
parameter namely the ratio of the radii of these two spheres R s /Rq = t s /tq. It is also 
assumed that the photon injection takes place impulsively at t = 0, at energies ~ 1 MeV. 
For the treatment of this more detailed model we have used a Monte-Carlo code developed 
by one of us (XMH). The details of this Monte-Carlo code have been discussed elsewhere 
and will not be repeated here. The interested reader is referred to the relevant publications 
(Hua & Titarchuk 1995). 

In Figure 3 we present the widths in time associated with photons escaping frorm the 
scattering cloud at several energy bands, as obtained by our Monte-Carlo calculation; the 
energy bands were chosen in the same fashion as in the analysis of Fenimore et al. (1995). 
The widths of the pulses as a function of the energy are given for three values of the ratio 
r s /ro, namely 1.0, 0.5 and 0.1. The photon injection was assumed to be monoenergetic at 
E = 1 MeV and the total Thomson depth of the source was taken to be To = 10. One 
can see that the value t s /tq = 0.5 provides good fit to the relation found by Fenimore et 
al. (1995), indicated by the straight solid line. 

The resulting photon spectrum, in vF v units, accumulated over the entire burst is 
given in Figure 4, with the spectra escaping at specific times given by the individual curves. 
The times are spaced equally in logarithmic intervals. The time-accumulated spectrum is 
consistent with that of GRBs in that it peaks at roughly 1 MeV. Clearly, injection of 
photons at higher energies and/or at larger Thomson depths could produce spectra which 
would simulate the observed ones more accurately. Finally, in figure 5 we present the rates 
of escaping photons as a function of time, for photons of a given energy range. The energy 
ranges were chosen to be identical to those used by Fenimore et al. (1995) to derive the 
energy - duration relation. 

III. Conclusions, Discussion 

We have presented above two models which can account for the recently discovered 
correlation of the GRB duration (or the duration of subpulses within a given burst) as a 
function of the photon energy. Both models appear to be able to reproduce the observed 
power law dependence of At on E given in Fenimore et al. (1995), as well as the general 
form of the observed GRB spectra. Considering in greater detail the comparison of our 
model spectra with those of GRBs (Schaeffer et al. 1994), the model spectra are generally 
consistent with those observed in that their luminosity (i.e. their vF v distribution) peaks 
at ~ 10 20 Hz, while exhibiting a power law low-frequency dependence, in agreement with 
observations. 

The slope of this low-energy power law is very well determined for the synchrotron 
cooling models, being proportional to i? 4 / 3 , the functional form of the single electron syn- 
chrotron emissivity. Interestingly, there exist several GRBs (e.g. GRB 910503 and GRB 
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910601) whose low frequency spectra are in agreement with this interpretation. In par- 
ticular, the spectrum of GRB 910601 (Share et al. 1994), seems to be fitted very well 
with the single electron synchrotron emissivity over its entire energy range, suggesting the 
synchrotron origin of the radiation associated with this specific burst. However, most of 
the spectra in Schaeffer et al. (1994) require low-frequency power law fits of significantly 
different slopes, typically vF v oc E 1 ! 2 . Clearly, the simple interpretation of synchrotron 
emission by a uniform, optically thin, impulsively injected electron population, such as de- 
scribed above, would be inappropiate for this type of bursts. However, more complicated 
models which relax the above assumptions could produce spectra agreement with observa- 
tion. For example, continuous (rather than impulsive) injection of the electron population 
over time scales long compared to the loss time scale of the electrons with Lorentz factor 
7 m would lead to spectra in agreement with this form. 

On the other hand, the downscattering model could accommodate spectra of the form 
vF v oc i? 1 / 2 , since the precise slope of the low energy power law depends on the distribution 
of the high energy photon sources within the cold cloud considered. These slopes are 
consistent with those of the time averaged spectrum of the downscattering model given in 
figure 4. In further support of this point, one should bear in mind that the spectrum of 
figure 4 represents the result of a 5-function photon injection at the highest energy bin, 
hence the upward turn near ~ 1 MeV in figure 4. Injection of a broader type of photon 
spectra or larger value for tq would result in smoother spectra, conforming more to the 
observation. 

In assessing the utility of the above models, one could claim that both may be of 
relevance in modeling the spectro - temporal evolution of GRBs. The synchrotron cooling 
model fares a little better in reproducing the duration - energy correlation in that it 
requires fewer parameters to do so than the downscattering model. On the other hand, 
the downscattering model, because of its additional freedom, appears to be able to provide 
more detailed fits to a larger number of GRB spectra. Discriminating between these models 
requires a more detailed combination of spectral and temporal analysis of individual bursts 
rather than the use of their average properties. Nonetheless, the form of the observed 
correlations along with their present model interpretation clearly indicate that, whatever 
their precise mechanism, the available energy of GRBs is originally released into the highest 
energy particles and their evolution is determined by its subsequent cascade to the lower 
energy ones. 

Despite the general agreement of the computed spectra of both models to those ob- 
served and their overall agreement to the duration - photon energy relation it is hard, 
at the present level of analysis, to draw any further conclusions concerning the physical 
mechanism underlying the emission of radiation in GRBs. This is due to the absence of 
an independent estimate of the parameters which set the scales of the time and the energy 
of peak emission in the vF v spectra, in either model. In the non-thermal model, the time 
is set by the magnitude of the magnetic field, or rather the combination B 2 ^ (or by p ra dl 
if inverse Compton is the dominant loss mechansim, and the bulk motion Lorentz factor 
T if relativistic expansion is of relevance), while in the downscattering model it is set by 
the density of the electron cloud through which the 7-rays propagate, both of which are 
indeterminate. The energy at which the h>F v spectra peak could be used to provide some 
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additional information on the radiative processes associated with GRBs. However, such 
information is independent, at this level of analysis, from that associated with the GRB 
time scales. In the non-thermal models, the energy of peak emission is set by the combi- 
nation 7 m 5, while in the downscattering models it is an independent parameter, defined 
by the injection of primary photons. 

Figures 2 and 5 exhibit an additional aspect of the spectro - temporal correlations of 
the GRB time profiles which has been touched upon only briefly so far, namely that of the 
time lags in peak emission among the various energy channels. This issue, only a peripheral 
one in this note, could be further developed both theoretically and observationally and 
might provide additional information on the nature of the GRB emission process, or at 
least a sufficiently strong criterion to reject one or both of the above models. 

Figure 5, indicates that the model of impulsive high energy photon injection leads 
invariably to lags in the peak (highest intensity) emission between the various energy 
channels, in the sense that the lower energy emission always lags behind the high energy 
one. Due to the non-dimensional character of the time axis of this figure, the lags have to 
be compared to the dispersion (FWHM) of the appropriate energy photons. It is apparent 
in the figure 5 the lag is smaller than the dispersion (FWHM) of the burst, in qualitative 
agreement with the results of fig. 12 of Norris et al. (1996). While the same data indicate 
generally lags significantly smaller than the observed FWHM of the pulses, we believe that 
more detailed study is required to decide the quantitative agreement of the model with 
the data and possibly exclude this model on the basis of this issue. 

In the relativistic electron cooling model, on the other hand, the presence of such lags 
is not ubiquitous, even for an impulsive injection of the electron population. As shown in 
figure 2, lags develop only for frequencies emitted by electrons with Lorentz factors lower 
than 7 m (to) ? as discussed earlier. Indeed, as seen in figure 2 there is little lag in the peak 
emission between frequencies ^ 10 20,5 Hz, while there are progressively increasing lags 
between the peak emission of lower frequencies. Since the time axis of figure 2 is absolute 
(i.e. in seconds) one can compare these lags direclty to those discussed in Norris et al. 
(1996). In doing so, one should bear in mind that the lowest channel of the above reference 
corresponds to curve / of this figure so that no lags longer than 0.5 sec should be expected 
for the chosen values of the parameters, in rough agreement with Norris et al. (1996). 

The issue of lags in peak emission within both models can be further complicated 
by the relation of the injection to the loss time scales in cases in which the latter is not 
impulsive. A further indication of the complex nature of this issue is the fact that, as 
noted by Norris et al. (1996), the lags in peak emission depend not only the energy band 
of observation but also the asymmetry in the burst time profiles, an issue not discussed 
here. We prefer not to delve into these issues in this note but concentrate instead on 
the relation between the observed photon energy and the burst duration, since that is, in 
both our models, property of the Green's function of the associated problem and hence 
independent of the details of injection. Some of the aspects of their findings could be 
accounted by the models presented above. However, since the entire burst time profile 
seems to be of importance in this issue, we will defer addressing it to future more detailed 
models of the spectro - temporal evolution of GRBs. The systematics provided by Norris 
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et al (1996) indicate that this issue is an important one and in need of further theoretical 
exploration. 

Finally, we would like to point out that neither model provides an answer to why GRBs 
emit most of their luminosity in the 7-rays rather than at other wavelengths. The down- 
scattering model treats the injection of primary photons as an initial condition and it is 
hence beyond the scope of any further discussion. The synchrotron emission model actually 
does fare a little better in this respect, since 7-ray emission is the natural outcome of elec- 
tron acceleration in shocks: It can be easily shown that for shock acceleration of electrons 
limited by synchrotron losses, the maximum Lorentz factor is of order 7 2 ~ (e//fi)(l/crr) 
and the frequency of the emitted radiation E 1 ~ (1/8/) (c/r ) ~ 2 x 10 22 Hz ~ (24//) 
MeV (r is the classical electron radius and / is a measure of the electron free path in 
units of the Larmor radius) , which has the interesting property of being independent of the 
magnetic field B. If this is indeed the reason that GRBs emit mainly in the 7-ray band, 
it would not leave much room for the relativistic blast wave models (Meszaros and Rees 
1992), since the latter postulate an additional relativistic boost of the produced photons 
with bulk Lorentz factors V ~ 10 3 in order that they are detected in the 7-ray band. 

The models presented above, despite their ability to provide good fits to the energy 
- duration relation of GRBs, they cannot, at the present level of analysis, provide an 
immediate, decisive answer to the question of the nature of GRBs. Nonetheless, the 
energy - duration relation reported by Fenimore et al. (1995) which represents the tightest 
correlation associated with GRBs to date, warrants their use in more detailed spetro - 
temporal analyses of GRB time profiles, which may set the ground work for uncovering 
the nature of GRBs. 

We would like to acknowledge useful discussions with M. Baring, J. Norris, J. Con- 
topoulos as well as useful, meaningful correspondence with the referee of our paper. 
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Figure Captions 



Fig. 1. -The evolution of the spectrum with time for the cooling synchrotron emitting 
relativistic electron model. The times corresponding to each curve are: (a) 10~ 3 s, (b) 
10" 2 s, (c) 10" 1 s, (d) 1 s, (e) 10 s. The B-field is assumed to be 30 Gauss, the maximum 
and minimum electron Lorentz factors at injection were taken to be respectively 10 9 and 
10 6 respectively. 

Fig. 2. -The synchrotron emission at individual frequencies as a function of time, for 
the relativistic cooling electron model. The frequencies corresponding to the various curves 
are: (a) 10 24 Hz, (b) 10 23 Hz, (c) 10 22 Hz, (d) 10 21 Hz, (e) 10 20 Hz, (f) 10 19 Hz, (g) 10 18 
Hz, (h) 10 17 Hz. 

Fig. 3.-The pulse width - energy relation for the downscattering model for different 
values of the t s /t$ ratio. 

Fig. 4. -The time-accumulated spectrum, for the downscattering model described in 
the text. 

Fig. 5. -The time evolution of emission for various energy channels, for the downscat- 
tering model. 
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